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1.  Introduction 


In  this  pspsr  w«  shall  consider  s  multi-item  inventory  probleat  with 
unspecified  single-item  unit  costs,  lather  than  examining  a  single  cost 
function,  we  shall  deal  with  an  approach  which  Incorporates  aggregate 
objectives  and  constraints.  The  objectives  are  low  investment,  low  total 
costs  and  high  service-level.  The  constraints  are  storage  room  capacity 
and  available  workload  for  handling  the  orders.  These  capacities  might  be 
increased  or  decreased  in  certain  fixed  quantities;  such  changes  in  workload 
and  storage  room  incur  costs  which  are  independent  of  whether  or  not  these 
capacities  are  used  to  their  full  extent.  The  objectives  are  conflicting 
and  in  many  real  world  problams  there  seems  to  be  no  single  cost  function 
for  determining  an  optimal  decision.  This  is  due  to  the  fact  that  there  is 
no  decision  unity;  Instead,  there  are  different  departments,  different 
managers  end  differing  Interests  involved.  How  much  the  company  should 
invest  in  Inventory  and  what  service  level  should  be  required  cannot  be 
determined  simply  by.  evaluating  a  single  cost  function,  but  is  rather  a 
result  of  Intensive  discussions.  What  operations  research  can  offer  is  a 
description  of  the  relationship  between  investment  and  service-level;  i.e., 
for  a  certain  Investment  we  are  able  to  maximise  the  service-level  under 
the  constraints  of  storage  room  and  workload.  Furthermore  we  might  study 
the  effect  upon  an  alteration  of  the  constraints.  This  enables  the 
management  to  find  an  "appropriate"  decision  weighing  the  different  interests 
ascertaining  that  there  is  no  solution  with  lower  total  costs  and  higher 
service-level . 

Despite  the  enormous  number  of  papers  on  Inventory  models  there  are 
only  a  few  articles  concerning  this  important  problem.  Most  papers  start 


by  assuming  that  marginal  holding*  shortaga  and  ordsring  costs  ara  givan. 

But  marginal  ordering  costs  ara  difficult  to  measure.  Most  suggested 
approaches  for  determining  ordering  costs  in  the  accounting  literature 
result  in  average  rather  than  marginal  coats  [15].  In  practice  there  ia 
usually  s  certain  workload  available  which  might  be  increased  in  fixed 
quantities.  Assigning  average  order  coats  therefore  does  not  solve  this 
particular  inventory  problem.  Holding  costs  should  be  composed  not  only 
of  the  cost  for  capital  but  also  of  marginal  costs  taking  into  considera¬ 
tion  a  storage  room  restriction.  The  use  of  shortage  costs  in  inventory 
theory  has  not  been  adopted  by  moat  practitioners  [2]  since  there  is  no 
basis  for  their  measurement  in  accounting  methodology  [15]. 

Only  a  few  authors  daal  with  an  aggregata  inventory  problem  as  described 
above. 

Starr  and  Miller  [13]  have  considered  an  "optimal  policy  curve"  for 
deterministic  demand.  Schrady  sad  Cbos  [11]  considered  a  continuous  review 
inventory  system  with  constraints.  Gardner  and  Daanenbrlng  [3]  extended 
Starr's  and  Miller's  approach  by  considering  a  continuous  review  stochastic 
model.  They  presented  a  method  that  avoids  cost  measurement  problems  and 
incorporates  aggregate  objectives  and  constraints.  They  describe  a 
procedure  for  obtaining  an  optimal  policy  surface,  the  axes  of  which  are 
measured  in  aggregate  terms:  the  percentage  of  inventory  shortages,  as  a 
measure  of  customer  service;  the  workload  in  terms  of  the  number  of  annual 
stock  replenishment  orders;  and  total  investment,  the  sum  of  cycle  and 
safety  stocks.  Aggregate  Inventory  decisions  are  defined  as  the  selection 
of  a  combination  of  the  three  variables.  This  model  obviously  reflects 
the  true  decision  problem  in  practice  much  better  than  a  single-item  cost 
model.  In  fact  data  simulation  is  frequently  employed  in  practice  to 
solve  decision  problems  described  above. 


Unfortunately  the  underlying  inventory  model,  considered  by  Gardner 
end  Dennenbrlng,  has  some  disadvantages  which  makes  It  inapplicable  in  many 
situations.  Host  of  the  Inventory  systems  Installed  in  practice  have 
periodic  review  and  not  continuous  review  [5].  It  hes  been  shown  that, 
given  a  certain  inventory  policy,  the  service-level  turns  out  to  be  very 
different  in  a  periodic  and  in  a  continuous  review  system  [91.  Further¬ 
more  we  should  be  aware  that  there  are  different  definitions  of  service- 
levels  [10],  which  result  in  very  different  ordering  policies. 

In  the  present  paper  a  similar  approach  to  that  of  Gardner  and 
Dannenbrlng  is  used.  But  there  are  some  essential  differences.  First,  a 
periodic  review  multi-item  inventory  model  life  considered.  Second,  we 
consider  only  two  objectives,  investment  and  service- level,  subject  to 
workload  and  storage  room  restrictions.  It  seems  to  be  more  realistic  to 
begin  with  a  given  workload  and  storage  room,  which  can  be  expanded  in  fixed 
quantities,  rather  than  to  assume  that  workload  and  storage  room  are 
continuous  variables.  Furthermore  an  overall  cost  evaluation  is  considered 
including  costs  of  investment,  storage  room  and  workload. 

An  interactive  algorithm  is  presented  which  allows  the  selection  of  a 
combination  of  service-level,  investment,  workload  and  storage  room  which 
is  "appropriate"  for  the  management.  This  method  produces  combinations 
which  lie  on  an  optimal  surface. 

Since  in  our  method  there  are  some  approximations  involved,  we  shall 


prove  the  validity  of  the  approximation  formulas  by  means  of  a  Monte  Carlo 


study  in  the  final  section. 
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2.  The  jjodel 


Consider  an  Inventory  model  with  n  items.  The  stock  of  every  item  is 
Inspected  et  the  beginning  of  e  review  period  and  an  order  is  placed  for 
those  items  for  which  the  stock  level  has  fallen  to  the  reorder  point. 

After  a  known  lead  time  X  the  orders  will  arrive.  The  demand  of  an  Item  in 
different  periods  is  a  random  variable  with  known  distribution.  Let  us 
define  for  item  k,  k-1,2,3 . . 


*kt 


-  inventory  on  hand  plus  on  order  at  the  beginning  of  period  t,  before 


an  order  is  placed. 


Q^t  -  order  which  is  placed  at  the  beginning* of  period  t 

-  stochastic  demand  in  period  t.  The  demand  in  successive  periods  is  a 

sequence  of  independent  and  identically  dlstrlbuged  random  variables 

2 

with  cdf  Ffc(r),  mean  end  variance  0^. 


pfc  -  price  of  item  k. 


Furthermore,  it  is  assumed  that  demand  which  cannot  be  immediately 
eetlsfled  is  backordered.  We  discuss  a  stationary  inventory  model  and  thus 
it  is  sufficient  to  consider  a  single  period  Inventory  model  where  the 
distribution  of  X^t  is  the  stationary  distribution  \|/(x)  16].  We  will  also 
assume  that  a  service  level  Y  •  which  is  defined  as 


Y 


1- 


cumulative  backorders  per  period 
average  demand  per  period 


is  an  appropriate  measure  of  customer  service  for  product  k«l,...,n. 
Notice  that  this  definition  of  a  service  level  is  equivalent  to  assigning 
shortage  costs  which  are 'dependent  on  the  amount  of  items  short  and  the 
length  of  time  the  shortage  lasts.  A  formal  proof  is  given  in  [9] .  This 
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type  of  shortage  cost  la  considered  by  aost  authors  [4],  [11],  [14]. 

In  whet  follows  we  consider  two  objectives: 

0^:  minimize  the  total  Investment  as  the  sum  of  cycle  and  safety  stock 

02:  maximise  the  service  level  y  and  two  constraints 
and  two  constraints 

C^:  the  number  of  stock  replenishment  orders  is  restricted  by  the  workload 
capacity 

Cji  there  is  a  storage  room  capacity  restriction 

This  two-objective-decision  problem  is  solved  by  determining  the  optimal 
policy  surface.  For  every  point  on  the  surface  it  holds  that  none  of  the 
objectives  can  be  improved  without  diminishing  the  other.  Before  proceeding 
we  shall  make  some  remarks  concerning  the  objectives  and  constraints.  The 
objectives  reflect  of  course  cost  considerations  and  the  constraints  are 
associated  with  costs  which  become  relevant  if  alterations  of  the  constraints 
are  allowed.  With  0^  we  only  control  the  variable  holding  costs  Induced  by 
Invested  capital.  If  the  storage  room  is  fixed  the  costs  for  holding  this 
room  are  fixed  too  and  thus  lrr el event  for  a  decision.  But  often  it  is  the 
case  that  storage  room  is  rented  in  certain  quantities  and  hence  the  costs 
for  holding  a  storage  capacity  becomes  relevsnt  for  a  decision.  The  same 
is  true  for  the  workload  restriction.  We  therefore  consider  a  second  set 
of  objectives  without  constraints. 

0^:  minimize  total  costs  involving  cost  for  workload,  cost  for  invested 
cspltal,  and  cost  for  holding  s  storsge  room. 

OjS  msxlmize  service  level  y 

Both  sets  of  objectives  should  bs  available  for  a  decision  process  for 
selecting  sn  " appropriate"  solution. 
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W«  shall  now  formally  define  our  objectives  and  eonatralnca.  Zt  is  well 
known  that  whan  given  a  single  Item  Inventory  model  of  the  type  described 
above  an  (stS)  policy  Is  optimal  [6].  Ha  will  thus  consider  (s^,S^) 
policies  In  our  approach.  Let 

ElI^ls^.S^]  -  expected  inventory  at  the  end  of  a  period  for  product  k 
E[QjJs^(S^]  -  expected  order  quantity  for  product  k 
E[HOk|VSkl  -  expected  number  of  orders  per  period  for  product  k 
E[Cl|{sk»Sk}]  -  expected  Invested  capital  In  Inventory 
E[MO|{sk>Sk}]  -  expected  number  of  orders  per  period  1 

E[SR|{sk,Sk)]  -  expected  storage  room 
Rote  that 

K[Cl|{sk.Sk}]  -  J  PkECZk|akiSkl 
k-1 

n 

*[SR|{sk.Sk>]  -  l  akE(Ik|sk,Sk} 
k-1 

where  ak  Is  the  unit  storage  room  for  product  k.  We  finally  formulate 
the  two-obj active-decision  problem  as 


min 

Uk.sk) 


{v*k} 


subject  to 


E[Sl|Uk,Sk)] 

BlY({»k»SK))3 


(1) 

(2) 
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where  (see  [9])  M(*)  end  a(*)  ere  solutions  of  the  renewal  equations 
z  z 

H(z)  -  F(z)+  /  M(z-t)dF.  <t)  m<z)  -  f(z)+  /e(z-t)f . (t)dt 
0  0 

and  L+  and  L  are  defined  by 

L+(x)  -  J  <t-x)f£+1(t)dt 
0 

/”(t-x)f^+1(t)dt  xi.0 

L“(x)  -  • 

At-x)f£+1(t)dt  x  <  0 
0 


Where  f£+1(t)  is  the  pdf  of  demand  in  lead  tine  plus  review  time.  Note 
that  -  S^-s^.  Although  an  exact  solution,  i.e.  the  determination  of  the 
optimal  policy  surface,  is  in  principle  possible,  we  will  not  recommend 
it  since  the  approximation,  provided  below,  will  give  such  excellent 
results  that  there  is  no  incentive  to  carry  out  the  extraordinary  high 
computations.  Our  approximation  is  based  on  empirical  results  [14]  that 
there  is  only  little  loss  of  optimality  if  the  optimization  is  separated 


8— 


In  determining  first  Dk  ■  s^~*|c  *n£*  sf tsrvsrds  s^.  Furthermore,  we  consider 
Che  expected  values  under  the  realistic  assumption  that  is  large  and  we 
are  thus  able  to  simplify  the  expressions  (5)  to  (7),  an  approach  which  was 
introduced  by  Roberts  [ 8  J .  Following  these  principles  we  derive  the 
simplified  expected  values 

D. 

R[lk|(sk,Sk)]  -  Dk(l-»i  ^  )+  sk  -  (Xj+D^+U-Y)^  (8) 

«*>kl<vV>  '  (9) 


’VV  *  1‘ 


/  (t-sk)2fk+1(t/dt 


2[V  u2k/2ukj 


(10) 


+  a: 


where  + 


'i 


The  expression  (8)  is  derived  in  the  appendix;  (9)  is  a  well-known  limiting 
theorem  of  the  renewal  function  M(D)  (12]  and  (10)  was  derived  in  [9]. 

Since  the  constraint  (4)  will  always  be  active  we  can  solve  this 
optimisation  problem  using  the  Lagrange  method.  Let 


L(Q^» ... »Qq»0) 


n 


n 


I  P.E[I.  |Q.  ]  +  P[  l  E[N0.  IQ.  1  -  WLR] 
k-1  *  *  K  k-1 


(11) 


A  straightforward  application  of  the  Lagrange  method  yields 
<*  -  lc«l . » 


(12) 


2 


(13) 


The  reorder  points  sk  can  now  be  calculated  in  a  second  step  for  various 
service  levels  y.  Ue  determine  sk  by  equation  (10). 

The  cycle  plus  safety  stock  for  product  k  and  fixed  service-level  is 

then 

BUJvV  •  V1'0*5  q*>  +  V(Xk+1)Uk+(1“Y)Uk  (14) 

''k 


and  thus  the  total  Investment  is 

E[ClM  -  ZPkE[lk|(sk,Sk)] 


and  the  expected  total  storage  room  is 
E[SR|y]  -  EajEllJU^)] 

which  are  both  functions  of  the  service  level  y.  The  latter  value  is  now 
compared  with  the  storage  room  capacity.  The  service  level  can  be  increased 
as  long  as  the  storage  room  constraint  is  not  active.  He  obtain  a  diagram 
which  shows  the  Investment  versus  service-level  y»  this  diagram  will  end  at 
Yf  which  marks  the  maximum  sarvlce-level  at  which  the  storage  room  constraint 
becomes  active.  In  addition  a  diagram  which  gives  the  total  costs  (for 
investment  plus  fixed  costs  for  the  storage  room  capacity  plus  fixed  costs 
for  the  workload  capacity)  versus  service-level  can  be  presented  (see 
figures  2  and  3).  In  contrast  to  the  unit  single  item  costs  these  fixed  costs 
are  easy  to  determine.  These  diagrams  are  the  basics  for  selecting  an 
"appropriate”  service-level  and  investment  by  the  management.  The  steps  for 
obtaining  the  diagrams  are  summarized  in  the  following  algorithm. 
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Yes 


Yes 


Figure  1:  Algorithm 


I  selection  of  nev  storage  room  capacity) 
1  no 

|  selection  of  new  workload  restriction 


no 


»j.  v.  y  y  y»j 


In  this  subsection  we  shell  consider  en  example  with  only  100  items 
to  describe  the  algorithm.  In  practice  the  service-level  of  items  varies 
from  product  to  product.  Some  products,  for  example,  are  more  vital  than 
others.  Usually,  however,  one  is  able  to  classify  the  items  and  put  them 
into  certain  groups  which  have  the  same  service-level  requirements. 

We  thus  define  m  groups  such  that 

Yj  ■  gjY»  J"1.2,3, . . . ,m 


let  n^  be  the  number  of  products  in  subgroup  j  then 


.  u  m  m 

Y  “  n  I  Yk  "  n  l  “jV 
k-1  *  n  J-l  2  2 


which  yields  to  the  restriction 


•  X  ■ 1 


Furthermore 


Y  <  l/««*  (g.)  (16) 

j  2 

The  problem  of  assigning  is  thus  reduced  to  determining  m  numbers  g^ 
which  fulfill  (15)  and  (16).  This  will  usually  be  possible  in  a  reasonable 


We  shall  now  consider  a  situation  in  which  the  management  has  to 
decide  how  much  to  invest  during  the  year  to  come  and  what  service-level 


should  be  required.  There  is  a  certain  workload  available  which  allows 
22  orders  per  week.  The  cost  of  this  workload  is  $115,000  per  year.  Free 
handling  capacity  cannot  be  used  for  other  purposes.  An  Increase  or 
decrease  of  workload  is  only  possible  in  quantities  which  are  equivalent 
to  4  orders  per  week.  The  cost  for  such  a  workload  quantum  is  $15,000 
per  year.  There  is  also  storage  room  available  which  allows  for  the 
storage  of  60,000  units.  The  cost  for  renting  the  storage  room  is 
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$130,000  per  year.  Increase  or  decreaee  is  only  possible  in  quantities  of 
10,000  units  whictt  will  cause  costs  of  $20,000  per  year.  The  following  data 
are  aesuaed  to  be  known 

seen  deaand:  ^  -  20f9.89(k-l)  k  -  1.2,..., 100 

2  2  2 

variance:  -  c  k  -  1,2,..., 100 

c  m  0.2,  0.4,  0.S,  0.8,  1.2,  1.8 

unit  storage  room 

of  product  k:  is  randomly  chosen  to  be  in  the  interval  [0.5, 1.5] 

lead  time:  is  randomly  chosen  to  be  {2, 3, 4, 5, 6, 7, 8} 

price:  Pk  -  500/p^  k  -  1,2 . 100 

service  level:  5  groups  of  20  products  with  the  service  levels  0.8y, 

0.9y,  y.  1-lY.  1.2y 

We  assume  a  102  Interest  for  the  Invested  capital.  The  algorithm  is  demon¬ 
strated  for  c  -  0.2.  The  management  assigns  an  interval  for  the  service- 
level  602  i  y  <  1002.  According  to  the  algorithm  the  following  diagram  is 
plotted.  Due  to  the  workload  restriction  for  handling  and  the  storage  room 
restriction  the  maximum  average  service  level  is  79.52,  thus  the  diagram 
ends  at  that  point. 


Figure  3:  Sei 


0 
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Figure  4:  The  effect  of  an  alteration  of  the  workload 

restriction.  1)  initial  situation,  II)  increase 
III)  decrease  of  workload  restriction. 


Service  level 


We  see  that  a  reduction  of  workload  (curve  III)  results  in  a  lower  available 
service-level  while  an  increase  allows  a  higher  service-level  than  79. SZ. 
This  is  caused  by  a  change  of  the  cycle  stock,  which  decreases  (curve  III) 
or  increases  (curve  II)  the  available  storage  room  for  safety  stock. 

Dlagraa  5  gives  the  total  costs  for  the  alternatives. 


Figure  6:  Service- level  versus  Investment  for  different  storage 
room  restrictions 


Service  Level 


Since  the  workload  restriction  Is  fixed,  the  cycle  stock  Is  constant  and 
only  the  safety  stock  Increases,  which  results  in  a  curve  ending  at  a 
service-level  of  82. 5%. 

Figure  7  gives  the  total  costs  of  the  three  alternatives. 


Altogether  we  have  9  alternatives;  the  maximum  service-levels  and  total 
costs  of  these  service  levels  are  presented  in  Table  1. 


Table  1:  Maximal  service- level  and  total  costs  for  various 
combinations  of  storage  room  and  workload  capacity 


Stone*  too* 
Capacity 

Hbriaoad'^, 

Capacity 

M,000 

Has.  Sort.  Total  Coat* 

60,000 

Has.  So  nr.  Total  Costa 

70,000 

Has.  Sanr.  Total  Coata 

IS 

tor 

216 

69.32 

237 

76.52 

mm 

22 

73.32 

79.52 

231 

S2.52 

272 

2# 

7*.  32 

243 

62,32 

266 

S32 

2S6 

19- 


The  four  minima 1  cose  combinations  ore  presented  in  Table  2;  figure  8  gives 
the  minimal  total  costs  versus  service-level.  Table  2  shows  that  for  our 
example  a  maximum  service-level  of  73. SZ  can  bo  reached  by  the  combination 
of  workload  (22)  and  storage  room  (50,000).  A  service-level  of  79. 2Z  can 
be  reached  by  increasing  the  workload  capacity  to  (26).  Xf  a  higher 
service-level  than  79. 5Z  is  required  we  increase  the  storage  room  capacity. 

Table  2:  Optimal  Combinations  of  Workload 
_ and  Storaae  loom _ 


Maximal  Obtainable 

Service-Level  Workload  Storage  Boom 


-73. 5Z 

22 

50,000 

-79. 5Z 

26 

50,000 

-82. 5Z 

26 

60.000 

-83Z 

26 

70,000 

The  total  costs  can  be  obtained  from  figure  8. 


Figure  8:  Service-level  versus  Total  Costs 
of  Optimal  Combinations 
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The  management  now  has  all  tha  nacassary  information  to  decide  what 
aervice-level  and  corresponding  investment  should  be  selected,  and  whether 
or  not  the  workload  and  the  storage  room  restrictions  should  be  altered. 

He  notice  that  with  the  selection  of  a  certain  point  on  one  of  the 
curves  presented  above  we  determined  the  optimal  inventory  policy  as  well, 
l.e.  the  set  {s^.S^h 

3.  Validation  of  the  Model  by  Monte  Carlo  Simulation 

The  purpose  of  the  simulation  provided  in  this  section  is  twofold. 
First,  we  want  to  use  the  simulation  results  to  test  the  accuracy  of  the 
approximation  formulas  derived  in  section  2^-_  Second,  we  investigate  the 
variance  of  the  various  measures  of  system  performance  such  as  handling 
used.  Inventory  on  hand.  Invested  capital  over  the  periods.  Since  the 
restrictions  are  met  by  expected  values,  we  have  to  be  careful  that  the 
variance  is  not  too  large. 

We  have  performed  6  simulations  for  different  demand  structures,  i.e. 
c  -  0.2,  0.4,  0.6,  0.8,  1.2,  1.6.  For  C£  0.5  a  normal  distribution  was 

e 

used;  otherwise  the  gamma  distribution  was  found  to  be  appropriate  [9] . 
Notice  that  for  increasing  c  the  demand  becomes  very  erratic.  The 
simulations  were  run  for  workload  WL  -  22  and  service  level  y  •  0.82.  The 
theoretic  results  for  inventory  on  hand  and  invested  capital  are  given 
in  Table  3. 


Table  3:  Theoretical  results  (y  -  0.82  required)  for  mean  storage  room 
and  invested  capital. 


-  o/y 

Expected 
Storage  Room 
Used 

Expected  Invested  Capital 
In  Inventory 

In  $ 

0.2 

67.312 

68,855 

0.4 

86,185 

87,167 

0.6 

120,282 

120,445 

0.8 

162,687 

163,212 

1.2 

279,741 

283,268 

1.6 

440,667 

451,172 

Table  4  gives  the  simulation  results  with  1000  periods  and  50  repetitions. 


Table  4:  Simulation  results  for  mean  storage  room, 
mean  invested  capital,  mean  handling  and 
mean  service- level. 


Stores*  too* 


e  -  O/u 


Handling 


6S.322 

60 

86.761 

47 

122.373 

61 

166,104 

187 

286,733 

226 

443,082 

in. 

XovMMd  Capital 

Mean  Std 

Sarrlca-Laval 

Moan 

Mm*  - 

0.2 

49.838 

31 

+ 1.6 

822 

0 

0.4 

87,830 

43 

40.8 

S2S 

0 

0.6 

122.833 

S3 

♦2.0 

82.32 

0.3 

0.8 

167,129 

113 

♦2.4 

842 

2 

1.2 

289.783 

183 

♦2.3 

83.32 

1.3 

1.4 

432.373 

222 

40.2 

862 

4 

AX  1*  pareaataga  of  torUtiN  free  thwnUe  win 


tV 


WlTffy  >  f  l  »Ji  i  J!  y.  tj»j  y,  ij.  MJlLJAW^fW  jj  f ,  yj  nj/. 
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The  mean  values  of  inventory  and  Invested  capital  are  about  2Z  higher  than 
the  expected  values.  The  mean  workload  meets  the  constraint  for  normal 
demand,  while  for  sporadic  demand  c  >  0.5  the  actual  mean  workload  tends 
to  be  lower  than  the  constraint.  The  service-level  y  is  as  expected  if  the 
demand  is  normal.  For  sporadic  demand  the  actual  service-level  is  higher 
than  the  theoretical  value.  Tables  5  and  6  give  the  service-levels  in 
the  S  groups  for  c  ■  0.2  and  c  ■  1.6,  respectively. 

Table  S:  Service-level  y  for  c  ■  0.2 


VI 


Group  N 

Required  yZ 

Mean  of  Simulation  yZ 

1 

65.6 

65.8 

2 

73.8 

73.8 

3 

82.0 

82.0 

4 

90.2 

90.2 

5 

98.4 

98.4 

Table 

6:  Service-level  y  for  c  ■  1.6 

Group  N 

Required  yZ 

Mean  of  Simulation  yZ 

1 

65.6 

.  73.4 

2 

73.8 

80.4 

3 

82.0 

86.8 

4 

90.2 

92.9 

5 

98.4 

98.7 

The  expected  value  as  a  measure  of  performance  might  not  be  satisfactory 
alone.  In  some  periods  the  actual  values  fall  below  the  restrictions; 
in  other  periods  above.  But  if  the  variance  is  not  too  high  this  is 
acceptable.  We  therefore  will  present  the  standard  deviation  of  the 
actual  Inventory  on  hand, Invested  capital  and  handling  in  a  single 
period. 


We  would  expect  that  as  demand  becomes  more  erratic  the  fluctuation 


of  the  mentioned  variables  increases  and  thus  results  in  a  higher 
standard  deviation.  But  as  we  see  from  Table  7,  while  the  standard 
deviation  of  demand  increases  from  0.2(1  to  1.6(i,  the  standard  deviation 
of  the  inventory  on  hand  in  a  single  period  increases  only  from  0.1  x 
inventory  on  hand  to  0.75  *  inventory  on  hand. 


Table  7:  The  effect  of  demand  structure  on  the  performance  of  storage 
room  and  invested  capital. 


Storage  Room  Used 

in  1000  Units  Invested  Capital  Workload 


•  a/u 

Mean 

Std 

Std/Mean 

Mean 

Std 

Std/Mean 

Mean 

Std 

Std/Mean 

0.2 

68 

7 

0.1 

70 

6 

0.08 

22 

4.0 

0.18 

0.4 

87 

8 

0.09 

88 

7 

0.08 

22 

4.0 

0.18 

0.6 

122 

11 

0.09 

123 

10 

0.08 

22 

4.1 

0.19 

0.8 

166 

17 

0.1 

167 

15 

0.09 

22 

4.1 

0.19 

1.2 

287 

43 

0.15 

290 

39 

0.13 

22 

4.3 

0.2 

1.6 

445 

64 

0.15 

452 

57 

0.13 

22 

4.4 

0.2 

To  obtain  a  more  precise  picture  of  the  actual  performance  of  the  inventory 
system  under  consideration  we  studied  the  frequencies  of  inventory  and 
handling.  Table  8  gives  the  frequency  of  inventory  on  hand  in  a  single 
period  measured  in  terms  of  deviation  from  the  theoretical  expected  value. 


Table  8:  Frequency  of  Inventory  on  hand 


oh 

80Z 

90Z 

100Z 

110Z 

120Z 

130Z 

140Z 

150Z 

0.2 

0.04 

0.21 

0.39 

0.28 

0.07 

0.01 

0 

0 

0.4 

0.04 

0.23 

0.42 

0.26 

0.05 

0 

0 

0 

0.6 

0.02 

0.18 

0.47 

0.29 

0.04 

0 

0 

0 

0.8 

0.07 

0.14 

0.51 

0.31 

0.03 

0.0 

0 

0 

1.2 

0.003 

0.14 

0.62 

0.27 

0.006 

0.002 

0.002 

0.01 

1.6 

0.001 

0.18 

0.71 

0.09 

0,004 

0.002 

0.002 

0.01 

i 

The  conclusions  we  draw  from  Table  8  are  not  what  we  expected.  For 
c  «  0.2,  normal  demand.  In  39Z  of  the  periods  the  actual  Inventory  Is  as 
predicted  (100Z  of  expected  value).  In  21Z  of  the  periods  the  Inventory 
Is  10Z  below  the  expected  value,  while  in  28Z  of  the  periods  the  inventory 
is  10Z  above  the  expected  value.  As  the  variance  of  the  demand  Increases, 
the  concentration  around  the  expected  values  becomes  stronger.  For  c  ■  1.6, 
which  Is  a  very  erratic  demand,  for  71Z  of  the  periods  the  actual  Inventory 
Is  as  predicted.  We  notice  that  the  modal  value  of  the  simulation  corre¬ 
sponds  with  the  theoretical  expected  value.  The  distribution  of  the 
inventory  turns  out  to  be  shewed  to  the  right.  When  demand  becomes  erratic 
it  happens  that  In  a  very  small  number  of  cases  the  actual  Inventory  Is 
50Z  above  the  expected  value.  The  outliers  cause  the  increasing  variance. 

Table  9  gives  the  frequencies  of  the  orders  per  period.  These 
frequencies  are  almost  independent  of  demand  structure  given  by  the  value 
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Table  9:  Frequencies  of  orders  per  period 


c-o/u 

60% 

70% 

80% 

90% 

100% 

110% 

120% 

130% 

140% 

150% 

0.2 

0.02 

0.06 

0.12 

0.17 

0.28 

0.16 

0.11 

0.05 

0.02 

0.01 

0.4 

0.02 

0.06 

0.12 

0.16 

0.28 

0.15 

0.11 

0.06 

0.03 

0.01 

0.6 

0.03 

0.06 

0.12 

0.17 

0.28 

0.15 

0.10 

0.05 

0.03 

0.01 

0.8 

0.03 

0.06 

0.12 

0.17 

0.28 

0.15 

0.10 

0.05 

0.03 

0.01 

1.2 

0.03 

0.06 

0.12 

0.17 

0.28 

0.15 

0.10 

0.05 

0.03 

0.01 

1.6 

0.03 

0.06 

0.11 

0.16 

0.28 

0.16 

0.11 

0.06 

0.02 

0.01 

Cumula¬ 

tive 

0.03 

0.09 

0.20 

0.36 

0.64 

8r80 

0.91 

0.97 

0.99 

1.00 

The  reason  for  che  unexpected  result*  nsaely  that  fluctuation  decreases  as 
demand  variation  Increases.  Is  due  to  the  correlation  of  the  investigated 
variables.  The  following  figures  show  the  autocorrelation  function  of 
inventory  and  the  number  of  orders  for  c  •  0.2  and  c  »  1.6  respectively. 

We  see  that  for  normal  demand  we  obtain  cycles  of  high  inventory  and  a 
high  number  of  orders.  This  is  due  to  the  relatively  high  deterministic 
part  of  normal  demand.  It  is  obvious  that  during  every  4th  and  5th  period 
a  high  number  of  orders  arrives.  If  demand  is  sporadic  there  is  no 
significant  correlation  between  the  number  of  orders  in  different  periods. 
The  inventory  is  highly  autocorrelsted  as  one  would  expect  for  sporadic 
demand  since  there  is  no  demand  in  most  periods.  But  we  also  notice  that 
no  cycles  appear  if  demand  is  sporadic. 

We  might  conclude  from  our  simulation  study  that  the  performance  of  the 
inventory  system  was  close  to  that  predicted  by  the  theoretical  model. 

The  modal  values  are  the  same  as  the  expected  values  of  the  model.  But 
since  the  distribution  of  the  variables  under  consideration  tend  to  be 
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skewed  to  the  right  the  average  of  actual  Inventory  and  invested  capital 
is  2Z  higher  than  the  expected  value.  The  service-levels  are  as 
required.  The  variation  of  the  variables  are  higher  when  demand  has  a 
high  deterministic  part.  l.e.  c  <  0.5.  In  this  case  we  notice  the 

appearance  of  cycles  in  total  Inventory  and  orders.  It  seems  that  there 
has  been  very  little  attention  given  to  this  problem  up  till  now  in 
Inventory  literature. 
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9a  Correlation  coefficient  SO  of 
average  Inventory  on  hand, 
c  -  0.2 


Fig.  9b  Correlation  coefficient  SO.  of 
average  Inventory  on  hand, 
c  -  1.6 
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Fig.  9c  Correlation  coefficient  BO  of 
handling,  c*  0.2 


Fig.  9d  Correlation  coefficient  R0  of 
handling,  c  •  1.6 
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APPENDIX 


We  will  show  that  tha  expected  Inventory  on  hand  has  the  asymptotic 


value 


E(l|s»S]  -*•  D(1  -  +  s-Ux+1  +  (1-Y)U  as  D  «*■  » 


(Al) 


Note  that  we  omit  the  index  k  for  convenience. 


Proof:  Let 


S  0  D+s 

/(S-x)<j»(x)X+1d*f  J  /  (Dfs-y-x)<Kx)X+1m(y)dxdy 
S[l|a,S)  -  - - 0  0 - 


1+M(D) 


(1)  It  is  easily  seen  that 


/(S-x)$(x)X+1dx  ♦  Dfs  -  y. , .  as  D  -► 


X+l 


(11)  Note  that  (see  Smith  [12]) 


1+M(D)  ♦*  +  -4 
y  2]T 


l  as  D  •*>  • 


2  2 
where  Uj  •  U  +  a 


(Hi)  first  notice  thst 


D  D+s-y 


/  /  (Dfs-y-x)<Kx)X+1m(y)dxdy  •  /(D+e-y.  -y)m(y)dy-/  /  (D+s-y-x) 


X+l 


x  <Kx)X+lm(y)dxdy 


0  D+s-y 


(A2) 


then  the  first  term  at  the  right  hand  side  of  (A2)  is 


?  '  D  ^2  D  ^2 

(D+s-yx+1)M(D)-/ym(y)dy  -  <D+s-yx+1>(^  +^-|  -  l)+  D(^  + -  1) 


D2  U2 

-  fr  +  D(  -  i)  M  D  ■ 

M  2y 


The  second  term  of  (A2)  is  asymptotically  given  (see  [9])  by 


I)  00  _ 

/  /  (04s-y-x)$*+1(x)m(y)dxdy  ■*  ■—  /(x-a)2$(x)dx  as  D  +  * 

0  Ofa-y  w  • 

00 

/(x-a)2*(x)  -  (l-Y)V  •  114*1(0)1  and  with  (i),  (ti)  and  (iii)  we  obtain 
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